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ABSTRACT 

We study the averaging problem from a point of view of variation of spatial volume V . 
We show that in the space of spherically symmetric dust solutions which are regular on the 
spatial manifold the variation 5V vanishes at the Friedmann-Lemaitre- Robertson- Walker 
(FLRW) solution in an appropriate sense, which supports the validity of the FLRW solution 
as the averaged solution. We also present the second variation , giving the leading effect 
of the deviation from the FLRW solution. 
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I. INTRODUCTION 



The standard cosmology is based on the assumption that our Universe is homogeneous 
and isotropic. However, the present our Universe is not homogeneous but has clumpy 
structures like stars, galaxies, clusters of galaxies, and superclusters. The recent observation 

of the highly isotropic cosmic microwave background radiation is usually regarded as 
an evidence of the assumption, but this is the case only up to the stage of decoupling. 
Nonetheless, one may want to think that the homogeneous and isotropic universe model, 
known as the Friedmann-Lemaitre-Robertson- Walker (FLRW) model, reflects the averaged 
nature of the true Universe over a scale larger than a supercluster. In particular, we usually 
expect that the global expansion of the Universe is well approximated by the FLRW model 
with the energy distribution given by the volume average. 

However, this "averaging hypothesis" is never justified in a trivial way [0], since Einstein's 
equation is highly nonlinear. In particular, averaging over a (spatial) volume does not 
commute with the time-evolution in general, so the averaged initial data can develop in 
time in a quite different way from the true data. One is therefore required to clarify the 
meaning of the averaging and establish the domain of applicability of it. 

One of the earliest work connected to this problem is due to Futamase P], who built 
a formalism which gives the back reaction of small scale inhomogeneities to the global 
expansion along the spirit of Isaacson [^] using the post-Newtonian expansion. However, 
it is difficult to justify the approximation used in this approach, and the basic equations 
are still hard to handle. Zalaletdinov proposed a covariant averaging formalism starting 
from some axioms, but this is so complicated that it seems very difficult to draw useful 
consequences. Up to now, in spite of efforts including the above, the averaging has not 
been well understood, because of its complexities and conceptual difficulties. 

In this paper we study the problem in a quite different way from the conventional ones. 
We focus on the dynamics of the scale factor a(r), which is, if the space is closed |]TU 



equivalent to the dynamics of the total volume V{t) in the sense V{t) oc a^(r). We examine 
to what extent the FLRW solution can be a good model judging from a point of view of 
the time-development of the total volume V{t). Note that when given a space of solutions 
spanned by some arbitrary functions of space, we can think of V{t) as a functional on it. 
By evaluating the variation of V{t) at the FLRW solution we may be able to obtain some 
information about the quality of the FLRW solution as an "averaged" model. We explicitly 
do this for the spherically symmetric dust case. As a result, we find that the FLRW solution 
is a critical point for V{t) in an appropriate sense, which gives a good support for the validity 
of the FLRW solution as an averaged model. We also evaluate the second variation of V{t), 
which gives the leading effect of the deviation from the FLRW solution. 

The exact solution for spatially spherically symmetric dust spacetimes is known as the 
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Lemaitre-Tolman-Bondi (LTB) solution. We add to this solution the assumption that the 
spatial manifold be 5*^. We give a description of this subclass of the LTB solution first. 



II. THE SPHERICALLY SYMMETRIC DUST SOLUTION ON 

The metric is written with the synchronous and comoving coordinates as: 

d,2 = dr^ - e^(-'^)di?2 - r\r, R){d9' + sin' Od^'), (1) 

where X{t,R) and r{T,R) are the functions to be determined from the Einstein equation. 
The general solution to this metric is well known as the Lemaitre-Tolman-Bondi (LTB) 



solution (see e.g. [jTT],|T2|), which possesses three arbitrary functions f{R), F{R), and to{R). 
(Our notation follows Ref. fill, except for the sign of f{R).) With these the function A 
is given by e'^ = r'^/(l — f{R)), where the dash stands for the derivative with respect to 
R. The function r is given in three separate forms, depending upon whether the arbitrary 
function f{R) is negative, positive, or zero, and each solution possesses the FLRW limit of 
negative {k = —1), positive {k = 1), and flat {k = 0) constant curvature, respectively. (Such 
a limit is achieved when r{T,R) separates as r{T,R) = $(r)\I'(i?) [|12|-) Since our model is 
spatially S^, the spatial manifold does admit a constant curvature limit and it should be 
positive. So, it is natural to choose the positive sign of f{R) [TBI, for which the the function 
r is given by 

F F 
r = —(1 - COST/), r - To{R) = ^jy^iv - sinr/). (2) 

The arbitrary function to{R) is called the function of "big bang time", since each Killing 
orbit {R =constant) degenerates at time r = to{R). 

Now we have seen that the solution is parametrized by three "arbitrary" functions, but 
not all solutions are suitable for the spatial manifold S^. The aim of the remaining part of 
this section is to describe the conditions imposed on the three functions to obtain "regular" 
solutions on S^. (By a regular solution we mean that it is regular on the spatial manifold 
during a finite interval of time from the big bang.) We, by describing them, show that there 
does exist sufficiently large set of the regular solutions in the space of all formal solutions. 
This will be needed to make the formal calculations in the next section realistic. 

Let us think of the three-sphere as a sum of two balls, ~ U D^. This is similar 
to the decomposition of two-sphere into two disks, 5*^ ^ U D^, achieved by cutting the 
two-sphere along the equator. We can understand the spherically symmetric three-sphere by 
thinking of each ball spherically symmetric. Since each such ball has a symmetry center, the 
spherically symmetric three-sphere has two symmetry centers, which are degenerate Killing 
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orbits, as well. We label these points as i? = and i? = vr, so the relevant region of the 
coordinate i? is / = [0, vr]. The first condition we should impose is therefore 



r = 0, at -R = and vr. (3) 

The function r should be positive except at the boundaries. Furthermore, we must impose 
some regularity conditions. An efficient way to see this is to calculate the scalar curvature TZ 
and the curvature scalar polynomial TZabcdT^"'^'^'^ , where TZabcd is the curvature tensor. Using 
our local solution, we find 7^ = e and TZabcdR-"^'"^ = l2{F/r^f - 8{F/r^)e + where 
e = F'/lr'r"^) is the energy density of dust. Since these scalars must be finite, we have 

^ < oo, i? G [0, tt] (4) 



and 

F' 



<oo, Re [0,7r]. (5) 



While condition @ is necessary to avoid the conical singularity possibly generated at R 
and TT, condition (|) is imposed to avoid the well-known shell-crossing singularity [|14 



Finally, we have to impose the following coordinate condition, which is necessary for the 
coordinates to span the spatial manifold well; 

e^ = ^>0, i?e[0,7r]. (6) 

Without this condition, we possibly have fictitious solutions. 

To find the solution to the boundary conditions, note that there exists the freedom of 
reparameterizations R — 7(-R). Using this freedom we can fix the leading power of r at the 
boundaries to the unity, i.e., we can make 

r oc < ) \ . (7) 

[n-R (R^n) ^ ^ 

Near the big bang singularity (?7 — > 0), from Eq.(0) we have r ~ [F/ (4/))?7^, and t~to{R) ~ 
(F/(12/3/2))r^3, so we find r ~ (9/4)i/3iri/3(^ _ to{R)Y/^. Taking the condition (g) into 
account, we have 

FiR)ocl^' 3 . (8) 

Let 0" be a real number, and suppose rj oc R'^ {R 0) for an arbitrary fixed r. From 
Eq.(^, we find that if a < 0, the function r would oscillate heavily between positive values 
and zero near the boundary. So, this case is unsuitable for the condition. On the other 
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hand, if o" > 0, the boundary condition (^) is satisfied and the behavior of f{R) near the 
boundaries can be determined as 

^ ^ [ (tt - R)"^ (R^tt) ^ ^ 

where ai,a2 > 2. The condition (H) is now trivially satisfied. 

Next, consider the regularity conditions (^) and (^. Note that the function r on / 
for a fixed r has at least one extremum (FIG.Q), since r > on the interior of I and 
r = at the boundaries. From the condition (^, we find that the function F' should 
vanish where r' vanishes (FIG.^. Since F{R) is independent of r, r' should also vanish 
on / independently from r. Otherwise, the regularity would break instantaneously. A 
straightforward calculation gives 



1 

r = 



cos?7 



(10) 



where Aij]) = (1 — cos?]) — sin 77(77 — sinr^), and Bij]) = —(1 — cosr]) + | sin 77(77 — sin 77). 
Since r' is a homogeneous linear combination of /', F', and Tq with distinct coefficients as 
functions of r, the only possible points on / for which r' vanish independently from r are 
the points where F', /' and Tq vanish simultaneously. In particular, r' and F' should have 
the same sign everywhere on I since otherwise the energy density e would become negative 
in some regions. We exclude such a case for physical reason. On the other hand, the sign 
of /' is rather arbitrary, except that f{R) should be maximal and take value 1 where F' 
vanishes (FIG.^). Taking value 1 is a consequence of the condition (||). The sign of Tq should 
be opposite to F' or zero (FIG.^. This is a consequence of the condition that r' have the 
same sign as F' at least for a finite interval of time from the big bang (77 = 0). In fact, the 
leading powers of A{ri), B{ri), and sin 77 in Eq. (p!0|) are, respectively, j^r]'^, ~-^V^: so 
the term proportional to Tq dominates near the big bang singularity if Tq does not vanish. 
It is clear that if Tq had the same sign as F', r' would have the opposite sign when 7/ ^ 0. 

Now, we have spelled out all the conditions for F{R) > 0, f{R) > 0, and tq{R). They 
are Eqs.(|^) and (^, that f{R) should take maximum value 1 where F' vanishes, and that 
the sign of Tq should be opposite to F' or zero. It is nice to keep in mind that the profile of 
F{R) decides that of r if the solution is regular. 



III. VARIATION OF V AND AVERAGING 



Given an inhomogeneous solution, in what way can we measure the resemblance between 
the solution and the FLRW solution? We note that the only dynamical content of the FLRW 
model is the scale factor a(r), which can be regarded as the (cubic root of the) total volume 
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V{t) if the spatial manifold is closed [0. Since V{t) is always well defined for any spatially 
closed inhomogeneous spacetime [|I5[ there is a large amount of theoretical advantage to 
utilize it. Note that V can be considered as a functional V^[0i(x), ■ ■ ■0„(x)] on a space of 
solutions, where 0i(x), ■ ■ ■0„(x) are the arbitrary functions of space which span the space 
of solutions. For example, in the spherically symmetric dust case, n = 3 and we can put 
(0i(x), 02(x), 03(x)) = {F (R) , f (R) , To{R)) . If the inhomogeneous solution is not too far 
from the FLRW solution, we can compare the two solutions by evaluating the variation 
6V{t) at the FLRW solution. If this function 6V{t) vanishes for all r, we may say that 
the two solutions are dynamically close to each other and the "averaged" solution of the 
inhomogeneous solution corresponds to the FLRW solution. 

Below, we apply the above idea to the spherically symmetric dust solution on 5*^. Our 
space P of solutions are assumed to be spanned by the regular solutions on S^, and the 
variations are taken in it. More specifically, P is a subspace of the larger space P* which is 
spanned by all possible configurations of three smooth functions {F{R), f{R),To{R)), and 
defined by subtracting from P* the configurations which correspond to irregular solutions. 
The space P contains the FLRW solution, since this solution is regular. 

The total volume for the metric (|l|) is 

V{t) = Att Te^/VMP, (11) 





where c*" = r'^/(l — f{R)) and r is given by Eq.(^. For later use, we also define the total 
energy E: 

E = f edV = 4n T -j=^dR, (12) 

Js'^ Jo y/l- f 

which is a conserved quantity. Moreover, for convenience we change the parameterization 
{F{R),f{R),To{R)) to {A{R)J{R),To{R)) defined by A{R) = F/{2f/^). The standard 
FLRW limit is achieved in this parameterization when 

{A{R),f{R),Tom = (ao,sin2i?,roc), (13) 

where ao is a positive constant parameter, and tqc is another constant parameter. Since 
we can always choose tqc = by the coordinate transformation r — > r + tqc, parameter tqc 
is redundant as far as the FLRW solution is concerned, but we will find that in our wider 
context it is useful not to fix this parameter. At the limit ([13D , V is given by 

Vo(r) = 271^00^ (1 - cos?7)^, r - tqc = ao{r] - sin?]). (14) 

We vary the volume V with respect to A{R), f{R), and tq{R), and evaluate at the FLRW 
solution (ITSD- The formal result is 
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SV{t) =A7ral{l- cos 7]) Air]) ^ (35 A + tan RS A') sin"^ RdR 

Jo 

+2nal{l - cos rif T (tan di? 

^0 

/•TT 

— 47raQ(l — COS?]) sin?7 / {36tq + tanR6Tn)sm^ RdR. (15) 

JO 

We may expect that the second term vanishes if we take into account the boundary conditions 
in the previous section, but we have to check the continuity of the function tan RSf at 
i? = 7r/2 to do this. The zeroth variation of f{R) is sin^ and the maximal value should 
always be 1 as we explained in the previous section. This means that at the neighborhood of 
R = 7i/2, the significant (or "dangerous") variation of f{R) should always be approximated 
by sin^{R + a), where a is a parameter. Differentiating this with respect to a and putting 
a = 0, we find that 6f is approximated by 2 sin cos i? da, so that tan R6f is approximated 
by 2sin^i?da, which is continuous at i? = tt/2. We can hence safely omit the term. (In 
contrast to this, we cannot apply "integrations by parts" to the first and the third terms 
due to the discontinuities of sin^ RtanR6A and the similar term for 6tq.) Thus we have 

SV{t) =Anal{l- cos 7]) Air]) r {35 A + tan RS A') sin^ RdR 

Jo 

/•TT 

— 47rag(l — COS?]) sin?7 / (3(5ro + tan_R5rQ) sin^ _Rd-R. (16) 

Jo 

Note that the function rj depends only on r when it concerns the FLRW solution, so the 
functions of rj in the above expressions have been factored out of the integrals. Because of 
this feature 6V{t) becomes constantly zero, if (and only if) the variation of A{R) and tq{R) 
are taken so that the integrals in Eq.(|T^ vanish. This means that the FLRW solution is a 
critical point for the functional V{t), if we restrict the directions of the variation in that 
way. 

The meaning of the particular directions becomes clear if we consider the variation of 
the total energy E, which coincides with the first term of the rhs of Eq.(|l6D up to prefactor: 

6E = Stt {36 A + tan R6A') sin^ R dR. (17) 
Jo 

(We have dropped the zero term 127rao Jq {tan R6 f)' dR.) One may notice that the second 
term of the rhs of Eq. ([T6|) is given similarly by the variation of 

Jo Vl - J 

which has been defined by replacing A{R) by to{R) (and multiplying the factor —1/2) in 
the definition of E. In fact, we obtain 

SC = -A-K r {35to + tan RSt') sin^ R dR, (19) 
Jo 
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so that we can write the formula ([I6D as 

SV{t) = al{l - COST]) (^A{r]) ^ + sinr^^C^ . (20) 

Here, note that the space of solutions P is foliated by the surface Pe,c of constant E 
and C, while the FLRW solution O is a two dimensional subset in P spanned by the two 
parameters and tqc- We write the FLRW solution with fixed values of and tqc as 
Oao,TOc ^ We can easily find that there is a unique FLRW solution Oao,TOc ^^^h Pe,c, 
and we can define a map Av : P O hj this correspondence: 

Av : Pe,c ^ 0,,,ro.. (21) 

The significance of Eg. (pOD is that the FLRW solution Oaa,T0c that best matches a given 
inhomogeneous solution p E P is given by Av{p), the FLRW solution with the same energy 
E and the same "C". That is, since in each surface Pe,c the FLRW solution Av{Pe,c) is 
the critical point for the volume V{t), all the inhomogeneous solutions which are sufficiently 
close to Av{Pe,c') virtually manifest the same dynamical evolutions of volume as that of 
Av{Pe,c)- This is our main result. We will call Av the averaging map. 

The key relation ( [20| ) seems very reasonable (though it is not trivial) if we note the 



following fact about the FLRW solution. That is, if we vary the two FLRW parameters oq 
and the total volume Vo(t) for the FLRW solution will suffer a certain change. We can 
easily estimate it by differentiating Vq{t) with respect to the two parameters. (See Eq.(p^.) 
The result is 

dVoir) = Qn^alil - cos 77) {2A{^)daQ - sin r^droc) . (22) 

Moreover, using 

= 127r2ao, Cq = -67rVoc (23) 

for the value of E and C at the FLRW limit (|13D, we obtain exactly the same formula as 
Eq.®: 

dyo(r) = al{\ - COST/) [A{riy^ + siuT^dCoj . (24) 

In this sense the functionals E and C are "inhomogeneous generalizations" of the FLRW 
parameters oq and tqc. Relations (|23| ) also imply 



Av{Pe,c) = 0^^^^.^^^. (25) 



Recall that the parameter tqc contained in the FLRW solution is merely a kind of gauge 
freedom, and therefore we could fix it like tqc = 0. This seems to mean that we can restrict 
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the space of solution P to the shce C = by gauge fixing. This is in fact the case. Note that 
the coordinate transformation r — r — c, where c is a constant, induces the transformation 
To{R) To(-R) + c. Hence the space of solutions P contains gauge freedom of this type. 
Since the above transformation for tq(R) shifts the value of C by Anc Jq ^-^-^ d-R, every 

surface of constant C corresponds to the same set of spacetime solutions, and this shows 
the claim. Note that the gauge freedom we consider is the freedom of choosing the origin of 
the time coordinate. The restriction C = therefore provides us a natural way of choosing 
the origin of the time coordinate for every distinct inhomogeneous solutions, especially with 
non- simultaneous big bangs, viewing from the comparisons of V{t). 

We should comment on a specialty of the FLRW solution. It was natural to naively expect 
that the FLRW solution describes the averaged spacetime of an inhomogeneous spacetime 
that is homogeneous and isotropic over a large scale. We have seen that this is in fact the 
case for the spherically symmetric case in the sense 6V{t) = at the FLRW solution. One 
might similarly expect that a spacetime that is smooth over a large scale but fiuctuates over 
small scales could be approximated by a smoothed solution, or one might expect that a 
homogeneous but anisotropic (i.e., Bianchi or Kantowski-Sachs-Nariai type [0) spacetime 



solution could be served as a good model for an inhomogeneous spacetime that is homo- 
geneous over a large scale. However, at least for the spherically symmetric case only the 
FLRW solution is the good averaged or smoothed spacetime model, since there does not 
exist a point where the function ri{T, R) becomes a function of only time r in P, except at 
the FLRW solution, i.e., 6V{t) = holds only at the FLRW solution. 

The accuracy of the FLRW solution as the averaged one can be estimated from the second 
variation That is, we first parameterize the arbitrary functions with one parameter e 
so that e = corresponds to the FLRW solution, and expand them in power of e: A{R; e) = 
ao + e6A{R) + (l/2)e^5^A(_R) H — -, and similarly for f{R; e) and tq{R; e). Then, we expand 
V[A{R; e), f{R; e), tq{R; e)] in accordance with these expansions; 

V = Vo + e6V+ ^e^^V + ■ ■ ■ . (26) 

If the first order term vanishes as in our case, the second order term ^e^S^V gives the 
leading term for the deviation from the FLRW solution, i.e., the accuracy is given by v{t) = 
{V - Vo)/Vo ~ le-'S^V/Vo. 



After a lengthy calculation and a similar consideration as for 6V, we obtain 

5 V = al{l- cos r]) i A{r]) ^- + sin r] j 

+127rao(g{v)J[SA,6A] + (1 + 2 cosr^) J[5ro, (5ro] 

-2(3 sin ri-ri{l + 2 cosri))J[SA, Sto]) , (27) 
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where and are the second variations of, respectively, the total energy E and the 
conserved quantity C: 

6^E = 247r^^ (^6^ A + ^ tani?^^^'^ sin^ RdR 

+2Arr r ( + tan R6A6f' + (l-- sin^ r) 6f6A'] dR, (28) 

JO \ cos^ R cos'^ R \ 3 ) ) 

= -1271 (^5Vo + ^tani?5V^j sin^i^di? 

-127r r f ^ + tani?5ro5/' + ^ (l - ^ sin^ i?) 5f5r'\ dR, (29) 

«^ \ COS XL COS -iL \ o / / 

and we have defined 

J[A = [ ((■)(*) + ^tani? ((■)'(*) + (■)(*)')) sin^i^di?, (30) 

and G{fi) = 7 — 8 cos rj + cos^ rj — Grj sin rj + 17^(1 + 2 cos 77). (If we were allowed to apply 
"integrations by parts" , the above formulae would become much simpler, but because of the 
same reason as in the calculation of SV, we could not do so. A formal application would 
cause a divergence of the result.) 

Since our variations are taken in the surface Pe,c of constant E and C the variations of 
any order of E and C vanish, in particular, 6'^E = and = 0. Hence we obtain 



5V(r) = 12iTao[g{T])J[SA,6A] + {1 + 2cosr])J[6To,6To\ 

-2(3 sin// - ri{l + 2cosri))J[5A, Stq]). (31) 



This depends only upon the first variations of the functions A{R) and Tq{R). 

We close this section by giving an explicit example. The three functions are f{R) = 
sin^R, A{R) = 1/3+ (l/900)(sin3i? - 12sin5i? - 9sin7i?), and to{R) = 0. (This example 
has been made by taking f{R) = sin^ R and putting 

E{R)= r f\x)fHxM^)dx. (32) 







We can check that all the regularity conditions for E{R) are satisfied if iplR) is a pos- 
itive function on [0, vr] such that it makes E{tt) = 0. In particular, if we choose 
^(R) = 1 + (l/10)((l/3) sin 5/? - sin 7/2), we have F{R) = sin^ R{2/3 + (1/450) (sin 3i? - 
12sin5-R — 9sin7-R)), and the A{R) presented above.) The corresponding FLRW solu- 
tion as the averaged one is given by ao = 1/3 and tqc = 0. FIGs.|| to |^ are, respec- 
tively, profiles of the metric components r and e^, and those of the energy density e. 
FIGM shows the time-development of the total volume V. In each figure, the profiles 
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of the corresponding FLRW solution are also depicted with dashed curves. FIG.|] shows 
logj^Q Y|5^y|/Vo, which is the estimation of the accuracy log^g l"^! by the second variation 
S^V , and the dashed curve shows the exact one log^g \y — Vo|/Vo. The second variation 
S^V is evaluated by putting ebA = (1/900) (sin 3i? — 12sin5i? — 9sin7-R), so the integral 
e^J[574, (JA] = — ll7r/324000 ~ —10"'^. We can see that the accuracy v is in this example 
better than 10~'^ throughout the expansion phase, though the energy fluctuation becomes 
larger than 10~^. 



IV. CONCLUSIONS 



We have seen that the FLRW solution is the critical point for the volume V{t) in the 
space of spherically symmetric dust solutions on S^. In accordance with the fact that the 
FLRW solution contains two parameters (, though one of which is redundant in a sense), we 
found that there is a natural foliation in the space of solutions defined by constant energy 
E and another quantity "C" . The exact statement of our result is that in each leaf of the 
foliation there exist a unique FLRW solution and this point is critical with respect to the 
variations taken in the same leaf. In our view, the "averaged" solution for all inhomogeneous 
solutions in a leaf corresponds to the FLRW solution in the same leaf. 

Although our discussions have relied on the known exact solution for the spherically 
symmetric case, we have seen that the correspondence between Eqs. (PPj) and ( P^ is very 
natural and seems to be independent of the spherical symmetry. In fact, we have already 
obtained a preliminary result that supports a direct generalization of the present result [jl8|. 
This will appear elsewhere. 



We could not discuss a relation to observables like the distance-redshift relation [|T9 
which will be worth investigating further. 
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FIG. 1. A possible profile of r. There are three extremal points R = Ra, Rb, Rc in this example. 




FIG. 2. A possible profile of F. The extremal points decide those for r. 




Ra Rb Rc r 

FIG. 3. A possible profile of /. At the extremal points of F, f should also be extremal, and at 
the maximal points among them / should take value 1. 



FIG. 4. A possible profile of tq. The sign of the derivative should be opposite to that of F or 
equal to zero. 
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1.2 



FIG. 9. Accuracy of Vq{t) as the averaged solution — evaluated from the second variation 
logio tI*^^^I/^^o solid curve and from the exact one log^g \V — VqI/Vq for the dashed curve. 
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